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We determine the conditions under which topological order survives a rapid quantum quench. 
Specifically, we consider the case where a quantum spin system is prepared in the ground state of 
the Toric Code Model and, after the quench, it evolves with a Hamiltonian that does not support 
topological order. We provide analytical results supported by numerical evidence for a variety 
of quench Hamiltonians. The robustness of topological order under non-equilibrium situations is 
tested by studying the topological entropy and a novel dynamical measure, which makes use of the 
similarity between partial density matrices obtained from different topological sectors. 



Quantum states with topological order (TO) belong 
to non-symmetry breaking phases of matter that cannot 
be characterized by local order parameters The ab- 
sence of any local structure in these states is ultimately 
responsible for their robustness against certain forms of 
local noise. For this reason topologically ordered quan- 
tum systems with ground state degeneracy and persistent 
energy gap, such as the Toric Code Model (TCM), have 
been proposed as robust quantum memories Q. More- 
over, systems with TO can support non-Abelian anyons 
as quasiparticle excitations, in which case universal quan- 
tum computation based on topological, and thus more 
robust, interactions is indeed possible Qj. 

Evaluating the persistence of topological order under 
adverse conditions is a crucial step towards future ex- 
perimental realizations of topological quantum hardware. 
The robustness of the TCM has recently been studied 
with respect to noise 14 and thermal equilibrium at a 
finite temperature :S, The quantum phase transi- 

tion from topological to magnetic order in the TCM has 
also been analyzed [1, [ij [hJ. We can say that a sys- 
tem possesses TO if it is gapped and it has a degenerate 
ground state whose degeneracy depends on the system 
topology. If topology is to be responsible for the degen- 
eracy, it is expected that different ground states should 
be indistinguishable locally. We show here that this is 
certainly the case for the TCM and use this property to 
characterize the presence of TO. In the absence of a lo- 
cal order parameter, establishing the presence of TO is 
a daunting task. An effective order parameter for TO is 
the topological entanglement entropy [HI, H3, St opo : if 
S'topo(IV')) = 0, the state \ip) does not have TO; while 
if St opo (\ip}) > 0, the state has TO. S t0 po has already 
given insights into the nature of the TCM and its transi- 
tions to non-topological phases @, 0, H, H[. By contrast, 
the dynamical properties of TO have remained largely 
unexplored. 

In this work we investigate the dynamical response of 
a state with TO to a sudden quench, a topic that has 
attracted a lot of attention recently both experimentally 



|13j | and theoretically [14(. If the dynamics of the sys- 
tem is governed by a Hamiltonian that does not support 
TO then it is crucial to determine the extent to which 
TO survives and to assess the nature of destructive per- 
turbations. In particular, we examine the case of a fast 
quench starting from the ground state of the TCM [15|. 
We determine the influence of such a quench on the TO 
by studying the dynamical behaviour of the block entan- 
glement entropy [16j | and the topological entanglement 
entropy, whose dynamics is calculated here for the first 
time. In addition we introduce a new dynamical measure 
of TO, which is based on the similarity between states 
obtained from different topological sectors. 

We start by considering a system of spin-i particles 
situated on the edges of a square lattice with periodic 
boundary conditions. The TCM Hamiltonian on such a 
lattice is Ho = - J2 S A * - J2 P b p> where A s = <8)/ es A/ 
is a star operator applied to spins on edges that meet 
on a vertex and B p = ®i<z p Zi is a plaquette operator 
on edges that form a square (X and Z denote the usual 
Pauli operators). These operators stabilize the ground 
state of the model, that is, |* ) = A s \^f ) = B p \^ ). 
There is fourfold degeneracy in the ground state corre- 
sponding to different topological sectors 0. This is be- 
cause the incontractible strings along the vertical and 
horizontal directions, denoted by w\ and w)2, commute 
with the Hamiltonian and map l^o) m to orthogonal 
sectors. The ground state manifold can thus be writ- 
ten as L = spanjw 



;|* >; i,j = 0,1} with |* ) = 

I C I 1 / 2 J2 g eG a \®)> where g is a string of star operators 
and |0) = 1 00 • • ■ 0) is a reference state (vacuum) in which 
all the spins are pointing in the +Z direction (equiva- 
lently, one may choose to define the vacuum in the X- 
basis and the strings g in the Z-basis). The Abelian 
group G of all star operators on an m x n lattice with 
N = 2mn spins is generated under tensor multiplication. 
Due to the periodic boundary conditions, \G\ = 2 mn ~ 1 . 

The lattice is prepared in the ground state l^o) of the 
TCM Hamiltonian Hq. At t = a different Hamiltonian 
H is suddenly switched on. The evolution is then given 
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by = exp(-iiJt)|* )- If K 10;)} is the spect: 

of iJ then ex.p(—iHt) = ^/=i exp(— ie/t)|0i)(^/| and 
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This is the most general form of the quenched state. 

Firstly we prove that there exist families of quench 
Hamiltonians that cannot disturb TO. To this end we 
consider that the eigenstates of H can be written as gen- 
eral spin flips x acting on the vacuum, i.e., \<j>\) = x|0). 
The spin flips define a group X under tensor multiplica- 
tion, {x G X\x = Xi ® Xj ® • • • }. To apply the argument 
also to the group Z of all ^/-operators, we take the initial 
state to be in a general superposition of the four topolog- 
ical sectors. Then the quenched state of Eq. (QJ becomes 
\m) = \G\- l/2 Y, g& G, Xl z X ,i,3 a ij expHerf]^) with 
Fij = (0\xiiv\iV2g\0) . The only surviving terms are those 
such that x\ = w\w 3 2 g = g lJ and therefore we obtain 
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geG,ij=0,l 



The entanglement between a contractible region A of 
the lattice and its complement, B, is given by the entan- 
glement entropy Sa = — Tr(/?A lnp^), where pa is the re- 
duced state corresponding to region A. Using the replica 
trick [13, S A (t) = -linwd„Tr[p^)], where p A {t) = 
Tr b[\^ '(t)) r (i)|] is calculated for the general quenched 
state. Specifically for the state of Eq. ((2]), we have 
P(t) = (■ L J2 g , g eG,i jk l v -«,,>;,£■>>■' ig.ggU/O gg. 
where £ l3kl (g,g) = exp[— ie(g lJ )t] exp[ie(g kl )t] and 
we have used the property g = gg of the group G. 
In order to calculate the state p A {t) = Trs[p(t)] we 
distinguish the region, |0) = |0a) ® |0_b), and we also 
define the subgroups of G that act trivially on the 
subsystems, i.e., Ga = {g £ G\g — qa ® Is} and 
Gb = {g G G\g = 1a ® ffs}- Since A is contractible, 
there are always group elements hi,h>x G G such 
that ZiiWi G Gb for i = 1,2. In other words, it is 
always possible to deform the loop Wi so as to move 
it out of the region A. Therefore, we have pA{t) — 
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a iJ \ 2 £^(g,gg)g A \0A)(0A\gAgA- 



Notice that for the ground state of the TCM, pa does 
not depend on the coefficients a^, meaning that two 
orthogonal ground states are locally identical. 

The particular structure of the Hamiltonian H im- 
plies that e(g lJ ) = k 13 e{g) where k lJ is a constant; hence 
p\(t) = |G|-" E flJ €G A €G A n (*) 9i,a\Va){Ua\9i,a~9i,a 
Ql^AQl^A ' ' ' 9n, a\0a) (0A\g n ,Ag n ,A, where the 
U(t) = £(g 1 ,g 1 gi) £ (52,5252) • • ■ £(g n ,9n9n) includes the 
time dependence of the state. Each expectation value 
of the form (0a\9i,a9i,a9i+i,a\0a) is set equal to 1 and 
this imposes the condition gi+i,A — 9i,A9l,A and so 
9i9i9i+i G Gb- In this way we get p A (t) = \G B \ n \G\~ n 
T, gun eG,g n £G A n (*) 9l,A\0A)(0A\9n,A9n,A, where now 
II(t) = £(gi, g n g n ). Using the replica trick we find 



S A (t) =]n(\G\/\G A \\G B \)- Therefore the entanglement 
entropy of the quenched state \^(t)) of Eq. $2$ remains 
constant in time. Moreover this state obeys an area law 
for the entanglement, S = La — 1, where La is the length 
of the boundary of A and the universal finite correction 
signals the presence of TO [l2[ • 

For the topological entropy we use the expression due 
to Levin and Wen 11], S top o = -Sabcd + Sabc + 
Sacd — Sac where A,B,C,D are the sides of a 2 x 2 
square in the lattice each containing two spins Q. Evi- 
dently, Stopo is time-independent for the quenched state 
of Eq. (0) since the block entanglement entropy is con- 
stant for this state. This completes the proof of the first 
result, which can be stated as follows: Stopo (i) = I if the 
eigenstates of the quench Hamiltonian can be written as 
spin flips on the vacuum (in either basis) . For this reason 
we call such quenches non- destructive. Any Hamiltonian 
H(Z) or H(X) that depends on only one Pauli operator 
has such eigenstates. Physically relevant cases include 
Hi = -J2iLi h iSi and H 2 = - J2{i,j) JiSiSj, where 
S = X, Z and i, j are nearest neighbors. In the absence of 
disorder (i.e., for hi = h and Ji — J) we have K^olV'W)! 
= I C| 1 2~2 g eG C0S [E{g)t] , where E(g) are the eigenen- 
ergies in either case. Using the commutativity of terms 
in H12 it is easy to show that such quenches have a pe- 
riod of recurrence T — n/2h,Tt /2J independently of the 
lattice size. For this reason, we say that such a quench 
is trivial Moreover, the entanglement remains constant 
at all times even in the disordered case. Conversely, a 
quench will be non-trivial if the periodicity is lost in the 
large system limit and the recurrence time increases with 
the system size. In order to detect the periodicity of the 
system we will study the overlap | (^0 1 ^C^)) I between the 
initial state |^o) and the quenched state of Eq. Ji}. 

In order to detect how TO is preserved by a non-trivial 
quench we use St opo and, in addition, the Uhlmann fi- 
delity T{t) between reduced systems that correspond 
to orthogonal ground states. The T{t) is given by 

Hi) = Tr^i(<K(J)^W, where a* (t = 0) is 

the density matrix of u^u^l'I'o)- As we have seen above, 
in the topological phase the reduced systems are indistin- 
guishable and J-(0) — 1, by construction. In numerical 
calculations with finite size systems we need to choose 
the region A such that its bulk contains more degrees of 
freedom than the boundary. For this reason we use this 
measure only for maximum available lattice sizes. 

We are now ready to examine non-trivial quench 
Hamiltonians for small lattices by means of exact nu- 
merical diagonalization. The quantities of interest are 
I (*o I *(*))!, Stopo (t) and Tit). We start with 



J 2 2J XiXj. 



(3) 



For Ji 3> J2 this Hamiltonian reduces to Hi, while for 
J2 3> J\ it reduces to H^. We have shown that for both 
of these extremes the TO of the initial state is robust 
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FIG. 1: (Color online) Numerical results for the quench H3 
with couplings J1/J2 = 0.33. (a) Absolute overlap | ( M/ 1 ^ } j 
on lattices with N = 8, 12 spins, (b) Stopo (t) for the N = 8, 12 
lattices and fidelity J-{t) for the N = 12 lattice. 
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FIG. 2: (Color online) Numerical results for the quench H4 
with h = 0.34. (a) Absolute overlap |(* |*(t))| on lattices 
with N = 8, 12, 18 spins, (b) Stopo(*) for the N = 8, 12, 18 
lattices and fidelity J-{i) for the N = 18 lattice. 



during the evolution. By contrast, in the intermediate 
regime, the recurrence time increases with the size of the 
system. The two competing terms affect both loop struc- 
tures that support TO. The Stopo and the T{t) quickly 
drop to low values (we note that T(t) is typically more 
sensitive numerically than Stopo)- The results are pre- 
sented in Fig. [TJ The system undergoes an eigenstate 
thermalization in which TO is destroyed. Having used 
several values of J\ j J2 we find that the most disrupting 
case for TO occurs for J1/J2 ~ 0.33, which is the pre- 
cursor of the quantum critical point for H3. In Ising-likc 
models the speed of interactions is ~ \J J\Ji [II] and is 
maximum at the critical point. At this value correla- 
tions spread faster in the lattice and the system samples 
a greater number of states in the Hilbert space. 

Are there non-trivial quenches that preserve topologi- 
cal order? At first sight this would seem unlikely because 
in such cases the total state typically explores large parts 



FIG. 3: (Color online) Numerical results for the quench Hs. 

(a) |(* |*(i))i on the = 8, 12 lattices for Ji = J 2 = 0.033. 

(b) Stopo (t) for the N = 8, 12 lattices and Tit) for the N = 12 
lattice for Ji = J2 = 0.033. (c) The same as in (a) but for 
Ji — J2 = 3.3. (d) The same as in (b) but for J\ = J2 = 3.3 

of the Hilbert space. However the Hamiltonian 

s p i 

displays interesting behavior. A reasonable scenario is 
one where the system is prepared in the ground state of 
H4 for h = and then it is let to evolve for some h 7^ 0. 
This can happen, e.g., if the ground state of the TCM 
is precisely prepared using syndrome measurements but 
there is some defect in the laboratory implementation. 
The model was analyzed in Ref. [U and it was shown 
that there is a precursor of a quantum phase transition 
even for small systems (for N = 8 we have h c w 0.34). In 
particular, it was shown there that for h < h c the system 
is TO and for h > h c it is normally ordered. 

We find that this quench is non-trivial: as the system 
size increases the overlap | ('Jol^^)) I becomes smaller 
and recurrences occur later. In Fig. [2] we present the 
results for h = 0.34. The influence of this quench on the 
system is milder than H3 due to the fact that i/4 couples 
the initial state to a much smaller Hilbert space as it is 
constrained by B p = 1 for every plaquette p. On the 
other hand, the time-averaged Stopo increases with the 
system size. For small values of h this is in line with the 
statical behavior analyzed in Ref. Q: the ground state 
of i?4 for h 1 has TO and Stopo approaches the max- 
imum value in the thermodynamic limit. Nevertheless, 
there is an important difference with the statical case in 
that large values of h destroy TO (the ground state of H 4 
for h ^> 1 is a trivial separable state). For the sudden 
quench, by contrast, large values of h preserve TO as the 
size of the lattice increases. 

We also present a non-trivial quench where the differ- 
ent possible phases co-exist in the dynamical behavior, 

h 5 = -J2 a ° -Y, b p - J iJ2 z > - JzJ2 XiX J- ( 5 ) 
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For Ji > J 2 > 0.33 or J 2 > J\ > 0.34 the system un- 
dergoes a periodic dynamic as with the quenches Hi , H 2 ■ 
If J 2 <C Ji ~ 0.34 and Ji <C </ 2 ~ 0.33 the system has 
a non-trivial dynamic behavior where TO is preserved 
as with H4. Finally we have the most interesting case 
when J 1 ~ J 2 (see Fig. [3|). Here we have two differ- 
ent phases: if J\ ^> 0.33, the quench is destructive; but 
if J\ <C 0.33, TO is preserved, which has been verified 
by checking that the time-averaged topological entropy 
increases with lattice size. 

In summary, we investigated how a topologically or- 
dered ground state reacts to rapid quenches. TO was de- 
tected by the topological entropy and a new dynamical 
quantity, the fidelity between reduced states belonging 
to different topological sectors. Clearly TO is preserved 
if the quench Hamiltonian has a ground state belong- 
ing to a topologically ordered phase. Nevertheless, we 
showed that quenches in a different phase can still pre- 
serve TO. On the other hand, we also showed that there 
is a class of quench Hamiltonian (see Eq.Q), whose dy- 



namics leads to a suppression of TO. The destructive 
dynamics results from Hamiltonians that are both non- 
topologically ordered and have a finite propagation of in- 
teractions through the lattice. It is remarkable that once 
a state with TO has been initialized, it can be more ro- 
bust with respect to dynamical evolution, as opposed to 
statical perturbations. It would be interesting to study 
the evolution of the TCM in the presence of a suitable 
noise model under various quenches. This would be a first 
step in constructing a dynamical decoupling protocol [l9[ 
for the protection of the TCM from thermal noise. 
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